Distinguishing Spontaneous Quantum Hall States in Graphene Bilayers 
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Chirally stacked TV-layer graphene with A'^ > 2 is susceptible to a variety of distinct broken 
symmetry states in which each spin-valley flavor spontaneously transfers charge between layers. In 
mean-field theory the neutral bilayer ground state is a layer antiferromagnet (LAF) state that has 
opposite spin-polarizations in opposite layers. In this Letter we analyze how the LAF and other 
competing states are influenced by Zeeman fields that couple to spin and by interlayer electric fields 
that couple to layer pseudospin, and comment on the possibility of using response and edge state 
signatures to identify the character of the bilayer ground state experimentally. 

PACS numbers: 73.43.-f, 75.76.-fj, 73.21.-b, 71.10.-w 



Introduction. — Bilayer graphene [l|, |2| and its thicker 
iV-layer cousins, chirally (or ABC) stacked multilayers |3|- 
LZI, have attracted considerable theoretical [7H1J] and 
experimental [15|420| attention because of their suscep- 
tibility to broken symmetries that are accompanied by 
large momentum space Berry curvatures and different 
types of topological order. In a continuum model mean- 
field theory, the ground state is@, [2l| an Ising layer- 
pseudos pin ferromagnet in which each spin-valley fla- 
vor is|8l4lOl| layer polarized. The quasiparticle Hamilto- 
nian in these states develops mass gaps that change the 
character of the wavefunctions at small momentum and 
producejTl, [ill Berry curvature. The integral of Berry 
curvature over a suitably defined region of momentum 
space near a given valley is nearly exactly quantized at 
±27r. This property can be interpreted [22| as saying that 
each valley contributes zLe^/h to the Hall conductivity 
with a sign that reverses with valley index and with the 
sense of layer polarization. States with total Hall con- 
ductivity le^/h evolve |23l| smoothly into quantum Hall 
ferromagnets with v ^ I in the presence of a perpendic- 
ular magnetic field. 

When spin is ignored only two different types of states 
can be distinguished, ones in which the K and K' val- 
leys are layer polarized in the opposite sense produc- 
ing a quantum anomalous Hall (QAH) statejTI, luJ, l24| 
with broken time reversal (T) symmetry and orbital 
magnetization [7|, and ones in which the two valleys have 
the same sense of layer polarization producing an inver- 
sion (I) symmetry breaking quantum valley Hall (QVH) 
state[7|, lll[ with zero total Hall conductivity. When 
spin is included, there are three distinct states with 
no overall layer polarization as sumixiarized in Table ID 
i) a QAH state with Hall effect contributions of the 
same sign for opposit e sp ins, ii) a quantum spin Hall 
(QSH) stated El, [H-Ull with opposite QAH signs for 
opposite spins, and iii) a LAF state[7| that has QVH 
states with opposite layer polarization signs for oppo- 
site spins. Among these possibilities, lattice mean-field 
theory calculations [l2| suggest that inter- valley exchange 
weakly favors QVH states in the spinless case and LAF 
states in the spinful case. In this Letter we analyze how 



all three states respond to Zeeman coupling to their spin 
and to electric-field coupling to their layer pseudospin 
degrees-of-freedom. We find that the Zeeman field re- 
sponse distinguishes QAH states from QSH and LAF 
states. In the LAF, the Zeeman field induces a non- 
coUincar spin state in which the components of the spin- 
density perpendicular to the field arc opposite in op- 
posite layers, while those along the field direction grow 
smoothly with field strength and arc identical. The three 
states respond similarly to an electric field between the 
layers, which can induce first order transitions at which 
the total layer polarization jumps. 

There is already some suggestive experimental evi- 
dence for spontaneous quantum Hall states in graphene 
multilayers that is consistent with mass gaps A '^ 2 — 8 
meV in recent studies of suspended bilayers'l5-19| and 
trilayers|20| . Since the gaps are seen[16i, JS, 2M o^ily 
at temperatures well below A/fc^ they appear to be of 
many-body origin. Moreover, measurements of bilayers 
in a perpendicular magnetic field B appear to show that 
both 1/ = ±4 and v = Q quantum Hall states can persist 
to zero- magnetic field 17|, Il8| , implying that spontaneous 



quantum Hall states with total Hall conductivity quan- 
tum number / = 0, 4 can be stabilized by interactions at 
5 = 0. 

Continuum model mean-field theory. — In single-layer 
graphene the band dispersion remains linear over a broad 
range of energy surrounding the charge neutrality point. 
When graphenc's honeycomb layers arc chirally stacked 
only two sublattice sites, one located in the top layer 
and one in the bottom layer, are not connected to near- 
neighbors in other layers and are therefore relevant at low 

TABLE I: SummaryQ of spin-valley layer polarizations (i or 
h), broken symmetries, charge (C) and spin (S) Hall conduc- 
tivities {e^ /h units) and insulator types for the three distinct 
states with no overall layer polarization. 
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energies. Hopping between these sites, e.g. from top (t) 
A to bottom (6) B, becomes an A^-step process, leading to 
two remarkably flat bands with -iik^ dispersion and layer 
pseudospin chirality iV[J,0|- These unique band features 
are encoded in the low-energy k-pHamiltonian given be- 
low. Because of the flat bands and the large pseudospin 
chirality, interactions become dominant at low energies 
in few- layer [4l [D] chiral graphene. In mean- field theory 
inversion symmetry is broken [9| within each spin and val- 
ley, leading in a contact interaction model to the follow- 
ing Hamiltonian: 



H 



HF 

ha 
h-p 



XI ^fc"s [^0 + /iH + hp] Ckps' , (la) 

kaf3ss' 

ek [ cos{N(bk)af + sm{N(bk)af] 4s' , (lb) 

[FoAor'^ + K.Az<^]4s', (Ic) 

-[Vo + V,arcT^J]Ki, (Id) 



where efc = (wo^fc) /(~7i) is the band dispersion, 

Vb,z = {Vs ± ^d)/2 denotes the average (difference) 
of intralayer and interlayer interactions, and A^| = 
A~^ '^/^{cj^gsi Ckas) i must bc determined sclf-consistently. 
Aq.z is the density sum (difference) of the top and bottom 
layers, cot (j)k = Tzk^/ky and Tz(±1) labels valleys K and 
K' . The Pauli matrices cr act on the which-layer pseu- 
dospin and s(±l) denotes the real spin. Because of the 
in-plane rotational symmetry of the continuum model, it 
is easy to verify that this mean-field Hamiltonian does 
not generate Hartree (H) or Fock (F) potentials that are 
off-diagonal in layer index. 

We seek self-consistent solutions for the N = 2 QAH, 
QSH, and LAF states. When Zecman coupling is ne- 
glected the Hartree and Fock contributions to the Hamil- 
tonian are mass terms proportional to a^, the four flavors 
decouple, and the mean-field equations are readily solved. 
For LAF, QSH, and QAH states the mass terms have 
the respective forms —msz ® cr^, —rriTz <Si Sz 'Si CTz, and 
—rriTz ® CTz where Sz is a spin Pauli matrix, as summa- 
rized in Table m Using the constant density-of-states per 
fiavor vq = 7i/(47r?i^UQ) of the normal state, introducing 
an ultraviolet cutoff at the inter-layer hopping energy 71 , 
and assuming weak-coupling, the gap equation can be 
solved to yield 



m = 271 cxp {-2/voVs) . 



(2) 



Influence of Zeenian Field. — When Zeeman coupling is 
included, the QAH state quasiparticles simply spin-split, 
leaving the ground state unchanged but the charge gap 
reduced. For a 4 meV spontaneous gap at zero-field, cor- 
responding to dimensionlcss interaction vqVs ^ 0.334 - 
close to the value expected to be appropriate for screened 
Coulomb interactions, a field of ~ 35 T drives the gap 
to zero. The QSH and LAF states, on the other hand, 
have more interesting non-coUinear magnetic-field in- 
duced states. We apply a Zeeman field in the x direction 
and allow spin-densities in the x — z plane. In practice 




Non-CoUinear LAF Order 

FIG. 1: (Color online) Upper panel: (Left axis) LAF tilt angle 
9 (green) and total effective-field tilt angle including both 
exchange and external field components (cyan) vs. in-plane 
magnetic field. (Right axis) Field aligned exchange mcosO 
(green) and total effective-field M+m cos (cyan) vs. in-plane 
field. Lower panel: sketch of the LAF tilt angles obtained 
from the upper panel. We assume a 4 meV spontaneous gap at 
5 = throughout the paper, corresponding to voVs ~ 0.334. 

this amounts to keeping A|^| real but allowing spin off- 
diagonal terms. In this case we find that for a 4 meV 
spontaneous gap, the LAF tilt angle 6 relative to the x 
direction decreases from 7r/2 at zero field to 7r/3 at 100 
T. The mass terms are correspondingly spin-dependent 
with components in the x and z directions. For the LAF 



h^^ = ho — m sin t 



[M + m cos 6] SxiEi ctq , (3) 



where 2M = g^L^B denotes the Zeeman splitting and m 
and 6 are determined by solving 



m sm f 



m sm V = 



mcosO 



AA ^ Es 

k.s=± " 



Vs V^ ^J^ + ITL COS 

n 2^ 



SCfc 



AA 



k.s=± 



E, 



(4) 
(5) 



with E±^^{M + m cos 6 ± euY + -m? sin^ 9. The four 
quasiparticle energies are :i^E±, so the gap is 2£'_ evalu- 
ated at tk = M + mcos9, i.e., 2msm9. 

For weak fields the quasiparticle spins are nearly per- 
pendicular to the Zeeman field. As the field strength is in- 
creased the quasiparticle state spin-polarizations, which 
are s and fc-dependent, all rotate toward the x direc- 
tion and the exchange field follows suit. Assuming that 
71 ^ m,M we find that the perpendicular LAF mass 




FIG. 2: (Color online) Canted LAF state quasiparticle bands 
for a series of in-plane magnetic field strengths: 10 T (blue), 50 
T (red) and 100 T (green). The dashed purple curve for 100 T 
shows the quasiparticle bands when the unaligned exchange 
field msin6' is neglected. 

component m sin 9 is still given by the right hand side of 
Eq. dH, and that 



mcos0 



i'oVs 



2 - lyoVs 
implying that the LAF tilt angle is 



-M. 



9 = arctan 



471 ■ (2 ^ ^0^3) 
gUBB ■ i^oVs 



,-2/(^oVs) 



(6) 



(7) 



This solution was confirmed numerically and is summa- 
rized in Fig [TJ 

The gap is nearly independent of M, in clear contrast 
to the QAH case. As Af increases the k = quasi- 
particle band extrema of the LAF move to larger k ex 
y/M + m cos 9 as illustrated in FiglJl For M ^ msm9 
the non-collinear LAF state can be viewed as an exci- 
ton condensate formed by pairing electrons in the bilayer 
majority spin band with holes in the minority spin band, 
fn this limit the LAF state is therefore similar to the 
Zeeman-coupling induced exciton condensate considered 
previously in the single-layer graphene case by Aleiner et 
al.^. 

Influence of Electric Field. — Because they all have a^ 
layer pseudospin order, LAF, QAH, and QSH states re- 
spond similarly to an electric field perpendicular to the 
layers, which adds a moCz term to the single-particle 
Hamiltonian. For the LAF, for example, the LAF masses 
m-f-4 for mo = differ only by a sign. When a perpen- 
dicular electric field is applied, masses are enhanced for 
one spin and suppressed for the other. In our mean-field 
calculations first order phase transitions occur between 
states with distinct broken symmetries as illustrated in 
Fig [31 leading eventually to a state in which the sense of 
layer polarization is the same for all spin- valleys |30[. Ex- 
perimental behavior in an external electric field will likely 
be sensitive to the pinning energies of domain walls that 
separate different spontaneous quantum Hall states. 



When an in-plane magnetic field and a perpendicular 
electric field are both present, the field aligned LAF or- 
der parameter m cos 9 is little changed compared to the 
E ~ case. The electric field dependence of rrif^i is 
mainly determined by a competition between moCTz and 
771 sin 6* Sz ^ Cz. The noncollinear LAF phase is, however, 
strengthened by its field-aligned order-parameter com- 
ponent and is more robust against a perpendicular elec- 
tric field when the Zeeman field is large, as illustrated in 
Fig [HI Fig [3| also shows that the LAF state stability can 
be dependent on the order in which the two fields are 
applied. Wc note that a small electric field between the 
layers can stabilize a state in which one flavor is polarized 
in a sense opposite to the other three and charge, val- 
ley, and spin Hall conductivities are all non-zero [7|. This 
state is not represented in Fig. [3] where we have assumed 
that the two valleys have the same layer polarization. 

Discussion. — Low-energy electrons in bilayer graphene 
have spin, valley, and layer two-component quantum de- 
grees of freedom. Because it appears in the band Hamil- 
tonian, the layer pseudospin plays a different role in bi- 
layer graphene physics than spin or valley. Flat conduc- 
tion and valence bands and Bloch states with J = 2 layer- 
pseudospin chirality combine to make the band state 
unstable toward a family of insulating broken symme- 
try states that have independent spontaneous layer po- 
larizations in each spin- valley component. Three dis- 
tinct states have no overall layer polarization, a quantum 
anomalous Hall state, a quantum spin Hall state, and a 
layer antifcrromagnet state [7|. In this Letter we have 
shown that the QAH state can be distinguished from the 
QSH and LAF states by examining the dependence of 
the charged quasiparticle gap on the strength of Zeeman 
coupling to an in-plane magnetic field. In the QAH case. 




Perpendicular Electric-Field Mass m [ meV ] 



FIG. 3: (Color online) LAF exchange-interaction masses rrif^i 
at parallel magnetic field B = T (magenta), 20 T (green) 
and 40 T (cyan) vs. perpendicular electric- field mass mo. 
The solid and dashed curves were obtained by following the 
state evolution vs. electric field at fixed Zeeman field and vs. 
Zeeman field at fixed electric field. The inner panel indicates 
the LAF (magenta) and fully layer polarized (green) state 
stability ranges vs. electric-field mass at zero magnetic field. 



the ground state is unchanged but the quasiparticle gap is 
reduced - vanishing when the Zeeman couphng strength 
is equal to the ground state gap via a mechanism rem- 
iniscent of the Clogston hmit in superconductors. The 
QSH and LAF states respond to Zeeman fields in a more 
interesting way, by establishing non-collinear spin states 
within each valley and evolving toward an unusual kind of 
exciton condensate in the strong Zeeman coupling limit. 
The gap of QSH and LAF states is independent of Zee- 
man coupling strength drawing a sharp distinction with 
the QAH case. When combined with probes that are 
sensitive to edge state transport, which is topologically 
protected[25[ in QAH and QSH cases but not in the LAF 
case, this property should enable any of the three states 
to be uniquely identified. 

It appears clear that bilayer graphene is exhibiting 
new many-body physics. This Letter points out that 
experimental studies of the Zeeman energy dependence 
of the gap could help to distinguish between different 
possibilities in bilayers, and also in larger N chiral few- 
layer graphene. As mentioned previously [2l| some the- 
oretical authors have concluded [13|, |l4| that the ground 
state of a neutral bilayer should be a nematic XY-planc 
layer-pseudospin ferromagnet which breaks in-plane ro- 
tational symmetry, rather than a 2-direction Ising pseu- 
dospin ferromagnet. (The z-component of the layer pseu- 
dospin density is the difference in density between the 
top and bottom layers while an x- or y-component indi- 
cates interlayer coherence.) The nematic states are most 
strongly distinguished from the X -symmetry breaking 
spontaneous quantum Hall states [7Hl2| by the absence 
of a charged quasiparticle gap in the former case. In the 
nematic state interactions generate mean fields that are 
off-diagonal in layer index and reduce the symmetry of 
the bands, splitting the 2tt K{K') Dirac points into two 
TT-Dirac points that are displaced from K{K') in an ar- 
bitrary direction. The mean-field-theory property that 
lower energy states are obtained with Ising compared to 
XY pseudospin order is related to the larger suscepti- 
bility associated with this pseudospin component. (The 
band eigenstates are perpendicular to the z-direction for 
all k, so all band states are easily rotated toward z pseu- 
dospin polarization.) Other potential explanations for 
the anomalies observed to date can be sought in trigonal 
warping effects, which are relevant below ^ 1 meV in bi- 
layers and have been ignored for simplicity in the present 
discussion, and in structural changes unintentionally in- 
duced by current annealing of suspended samples. There 
is however not yet a coherent explanation of how either 
of these might result in a gap at Dirac point. The ob- 
served gaps appears to be of many-body origin, in any 
event, since they appear only at t emp eratures that are 
much lower than observed gaps[16|-|20|. 
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